REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


p'jDiic  recx5rt»na  burden  ■^or  mis  ccHection  of  information  is  estimateo  to  averace  •  nour  oer  restxjnse,  mctuding  the  time  for  revtewmq  instructions,  searching  existing  data  sources, 
''atherma  ana  maintaining  the  data  needed,  and  comoieting  anc  reviewing  the  cci'ectton  or  mtormation.  Send  comments  regarding  this  burden  estimate  or  anv  other  aspect  of  this 
collection  of  mformation/mciuding  suggestions  for  reducing  this  burden,  to  Washington  Heaaauarters  Services.  Directorate  for  mtormaiion  Ooerations  and  Reports,  1215  Jefferson 
Davis  Highway.  Suite  1204,  Arlington.  VA  22202*4302.  and  to  the  Office  of  Managem'ent  and  Budget,  Paperwork  Reduaion  Projea  (0704-0188),  Washington.  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


REPORT  DATE 

9/5/96 


3.  REPORT  TYPE  AND  DATES  COVERED 

Interim;  4/10/96-9/30/96 


4.  TITLE  AND  SUBTITLE 

Moving  target  detection  and  motion  estimation  in 
foliage  using  along  track  monopulse  synthetic  aperture 
radar  imaain 


6.  AUTHOR(S) 

Mehrdad  Soumekh 


5.  FUNDING  NUMBERS 


N00014-96-1-0586 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

Research  Foundation  of  SUNY-Buffalo 
Office  of  Sponsored  Programs  Administration 
The  UB  Cormnons,  Suite  211,  520  Lee  Entrance 
Amherst,  NY  14228-2567 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORE5S(ES) 

Office  of  Naval  Research  Regional  Office  Boston 
495  Summer  Street,  Room  103 
Boston,  MA  02210-2109 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release. 


12b.  DISTRIBUTION  CODE 


13,  ABSTRACT  (Maximum  200  words)  ^ 

A  system  model  for  the  effect  of  uncalibrated  monostattc  and  bi static  radars  via 
gain  and  phase  functions  which  vary  with  the  radars  aspect  angles  and  the  radar 
frequency,  and  the  resultant  degradation  in  monopulse  SAR  images  and  their  difference 
which  is  used  for  moving  target  detection,  is  developed.  A  signal  subspace 
processing  of  monopulse  SAR  images  and  a  signal  subspace  difference  image  for  moving 
target  detection,  which  is  not  sensitive  to  the  unknown  gain  and  phase  error 
functions  of  the  uncalibrated  radars,  are  developed.  Numerical  results  are  provided. 


14.  SUBJECT  TERMS 


Along-tarck  monopulse  SAR;  uncalibrated  monopulse  radars; 
signal  subspace  processing. 


rs.  NUMBER  OF  PAGES 

25 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std,  Z39-18 
298-102 


Progress  Report  for: 

Contract  N00014-96-1-0586 

Moving  Target  Detection  and  Motion  Estimation  in  Foliage  Using 
Along  Track  Monopulse  Synthetic  Aperture  Radar  Imaging 

Submitted  to: 


Office  of  Naval  Research 

Program  Officer  William  J.  Miceli,  ONR  3131 

Ballston  Centre  Tower  One 

800  North  Quincy  Street 

Arlington,  VA  22217-5660 

Administrative  Grants  Officer 

Office  of  Naval  Research  Regional  Office  Boston 

495  Summer  Street,  Room  103 

Boston,  MA  02210-2109 

Director,  Naval  Research  Laboratory 
Attn;  Code  2627 
4555  Overlook  Drive 
Washingon,  DC  20375-5326 

Defense  Technical  Information  Center 
8725  John  J.  Kingman  Road 
STE  0944 

Ft.  Belvoir,  VA  22060-6218 


Distribution  Statement: 


Approved  for  public  release. 


Project  Director: 


Mehrdad  Soumekh 

Department  of  Electrical  &  Computer  Engineering 
State  University  of  New  York  at  Buffalo 
Amherst,  New  York  14260 
Phone:  (716)  645-2422,  extension  2138 
Fax:  (716)  645-3656 
Email:  msoum@eng.bufFalo.edu 
http:  / /www. acsu.bufFalo.edu/ ~msoum 


1 


pno  qTjALin  bssebos! 


Contents 


lA.  Overview  1 

iB.  Publications  1 

2.  Background  2 

3.  Effect  of  Uncalibrated  Radars  in  Monopulse  SAR  Imaging  2 

4.  Signal  Subspace  Processing  for  Uncalibrated  Radars  5 

4a.  Algorithm  5 

4b.  Implementation  and  Computational  Efficiency  7 
4c.  Stability  and  Noise  Effect  7 

5.  Future  Plans  9 

References  10 

Figures  11 

SF-298  25 


11 


Progress  Report: 

Moving  Target  Detection  and  Motion  Estimation  in  Foliage  Using 
Along  Track  Monopulse  Synthetic  Aperture  Radar  Imaging 


Mehrdad  Soumekh 

Department  of  Electrical  &:  Computer  Engineering 
State  University  of  New  York  at  Buffalo 
Amherst,  New  York  14260 


la.  Overview 

This  document  describes  the  progress  on  the  work  performed  for  “Moving  Target  De¬ 
tection  and  Motion  Estimation  in  Foliage  Using  Along  Track  Monopulse  Synthetic  Aper¬ 
ture  Radar  Imaging,”  under  Contract  N00014-96-1-0586  for  the  Office  of  Naval  Research 
during  the  first  six  months  of  this  contract  (4/10/96  -  9/30/96). 

The  scope  of  the  tasks  for  this  period  included  the  following: 

i.  Two-dimensional  system  and  signal  modeling  of  the  effect  of  uncalibrated  monostatic 
and  bistatic  radars  via  gain  and  phase  error  functions  which  vary  with  the  radars 
aspect  angles  and  the  radar  frequency,  and  the  resultant  degradation  in  monopulse 
SAR  images  and  their  difference  which  is  used  for  moving  target  detection; 

ii.  Developing  a  signal  subspace  processing  of  monopulse  SAR  images  and  a  signal 
subspace  difference  image  for  moving  target  detection  which  is  not  sensitive  to  the 
unknown  gain  and  phase  error  functions  of  the  uncalibrated  radars. 

This  report  includes  the  theoretical  derivation  of  the  gain  and  phase  error  functions  ef¬ 
fect  m  monopulse  SAR  images.  The  signal  subspace  modeling  and  processing  of  monopulse 
SAR  images,  and  constructing  the  signal  subspace  difference  statistic  are  described.  Nu¬ 
merical  results  are  provided. 


lb.  Publications 

The  work  has  produced  the  following  journal  article  which  would  contain  a  note  on 
the  support  from  the  Office  of  Naval  Research  under  Contract  N00014-96-1-0586: 

M.  Soumekh,  “Moving  target  detection  in  foliage  using  along  track  monopulse  syn¬ 
thetic  aperture  radar  imaging,”  accepted  for  publication  in  IEEE  Transactions  on 
Image  Processing. 
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2.  Background 

The  along  track  monopulse  SAR  imaging  system  utilizes  two  radars  for  its  data  col¬ 
lection.  One  radar,  Radar  1,  is  used  as  a  transmitter  as  well  as  a  monostatic  receiver.  The 
other  radar,  Radar  2,  is  used  only  as  a  bistatic  receiver.  In  our  original  proposal  and  [1], 
we  documented  a  signal  processing  algorithm  of  the  two  monostatic  and  bistatic  databases 
of  the  along  track  monopulse  SAR  system  to  obtain  two  coherently  identical  SAR  images 
of  the  stationary  targets  in  the  scene.  While  the  stationary  targets  appear  the  same  in  the 
monostatic  and  bistatic  SAR  images,  however,  the  same  is  not  true  for  moving  targets. 

This  fact  is  the  basis  for  developing  a  static,  which  we  refer  to  as  the  difference  image, 
for  moving  target  detection.  If  we  denote  the  monostatic  SAR  image  by  fmix,  y)  and  the 
bistatic  image  by  fb(x^y),  the  difference  image  for  moving  target  detection  is  defined  via 
the  following: 

fd{x,y)  =  fb{x,y)  -  fm{x,y).  (1) 

Figures  la-c  show  this  for  a  SAR  scene  in  which  the  two  radars  axe  fully  calibrated.  ^ 
This  example  corresponds  to  a  realistic  FOPEN  SAR  database  which  is  injected  with  the 
simulated  signatures  of  four  moving  targets.  The  parameters  of  this  SAR  scene  are  outlined 
in  our  original  proposal  and  [1].  Figures  la  and  lb,  respectively,  show  the  monostatic  and 
bistatic  SAR  images  of  the  target  area  which  contain  both  foliage  and  moving  targets. 
Figure  Ic  is  the  difference  image  which  shows  the  signatures  of  the  moving  targets;  the 
differencing  operation  has  removed  the  signature  of  the  stationary  targets  (foliage). 

3.  Effect  of  Uncalibrated  Radars  in  Monopulse  SAR  Imaging 

As  we  mentioned  earlier,  the  result  in  Figure  Ic  is  obtained  with  two  radars  which  are 
fully  calibrated;  i.e.,  there  is  no  relative  gain  and  phase  ambiguity  in  the  data  collected  by 
the  two  radars.  This  idealistic  scenario,  however,  is  never  encountered  in  practice.  In  a 
realistic  monopulse  SAR  system,  the  two  radars  exhibit  different  amplitude  patterns  (phase 
as  well  as  gain)  which  vary  with  the  radar  frequency  and  the  radar  position  (i.e.,  the  slow¬ 
time).  Moreover,  these  amplitude  patterns  vary  from  one  pulse  transmission  to  another 
due  to  heat  and  other  uncontrollable  natural  factors  which  affect  the  internal  circuitry  of 
the  two  radars.  These  subtle  changes  of  the  radars  amplitude  pattern  are  difficult  to  be 
detected  and  tracked,  and  are  unknown  to  the  user. 

To  develop  a  theoretical  model  for  the  undesirable  variations  of  the  amplitude  pattern 
of  uncalibrated  monopulse  radars  and  their  effect  in  the  difference  image  (for  moving 

^  Figures  are  attached  to  the  end  of  this  report. 
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target  detection),  we  begin  with  the  signal  model  for  the  monostatic  and  bistatic  SAR 
measurements.  We  denote  the  transmit-receive  amplitude  pattern  of  the  monostatic  radar 
in  the  spatial  (x^y)  domain  by  am{x,y  —  u,u))  where  (0, «)  is  the  radar  position  on  the 
synthetic  aperture  and  a;  is  the  radar  (fast-time)  frequency  [2].  Thus,  the  monostatic  SAR 
signal  for  a  stationary  unit  reflector  ^  at  (xn,yn)  is  [2] 

am{xn,yn  “  w,^)  exp[j2k^/xl  -F  {yn  -  m)^],  (2) 


where  k  =  u>/c  is  the  wavenumber. 

We  define  the  following  transformation  of  the  amplitude  pattern 

Atti  [2/?  sin  ^y7(a),  U?]  —  ^mij^Ti^yn  *^5^)) 

where 

^„(u)  =  arctan(^^^ - ),  (4) 

is  the  target’s  aspect  angle  when  the  radar  is  located  at  (0,  a).  (See  the  discussion  on  the 
Fourier  properties  of  the  S  AR’s  AM-PM  signal  in  [2,  Appendix]  for  a  physical  interpretation 
of  (3).)  Using  the  wavefront  reconstruction  algorithm  for  stripmap  SAR  [2],  it  can  be 
shown  that  the  Fourier  transform  of  the  monostatic  SAR  signal  of  (2)  with  respect  to  u 
(slow-time)  yields  the  following  target  function  in  the  spatial  frequency  domain: 

^m{kxiky')  —  Ajnikm^^  kxXn  jkyyn^j  (5) 

where  ku  is  the  spatial  frequency  domain  for  u,  and 

kx  =  -  kl 

ky  3=  k^l 

The  spatial  domain  monostatic  SAR  image  fm{^^y)  is  obtained  via  the  two-dimensional 
inverse  Fourier  transform  of  Fm{kx,ky). 

A  similar  analysis  can  be  performed  for  the  bistatic  measurement.  Suppose  the 
transmit-receive  amplitude  pattern  of  the  bistatic  scenario  is  ab{x,y  —  u,oj).  Then,  af¬ 
ter  the  phase  compensation  to  convert  bistatic  SAR  data  to  the  equivalent  monostatic 

^  One  could  consider  a  general  target  with  an  arbitrary  radar  cross  section  (amplitude 
pattern)  in  the  following  formulation  (see  [2]).  However,  for  notational  simplicity,  we 
absorb  this  amplitude  pattern  in  the  radar  amplitude  pattern. 
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SAR  data  ^  [1]  ,  one  obtains  the  following  SAR  signal: 


ah{xn,yn-u,u)  exp[j2k^/x^Tl^y^^^^^^]■  (7) 

Moreover,  the  bistatic  SAR  data  yields  the  following  spatial  frequency  reconstruction: 

Fb{kx,ky)  =  Ab{ku,uj)  exp(jkxXn+  jkyyn)^  (8) 

where 

Ab[2k  sin  6 niu),Lo]  =  ab{xn,yn  (9) 

and  (kx,  ky)  are  related  to  {ku,ix>)  via  (6).  The  spatial  domain  bistatic  SAR  image  fbix,y) 
is  obtained  via  the  two-dimensional  inverse  Fourier  transform  of  Fb{kx,ky). 

The  two  radars  are  considered  to  be  fully  calibrated  if 

ab{x,y  -  u,u!)  =  am(x,y  -  u,u}). 

In  this  case,  the  two  monopulse  SAR  images  of  stationary  targets  are  also  identical,  i.e., 

Fbi^kxy  ky^  —  Fffi{kxt  ky\ 

and,  thus,  the  difference  image  is  a  reliable  statistic  for  moving  target  detection.  This, 
however,  is  never  achieved  in  practice. 

In  a  practical  monopulse  SAR  system,  the  two  monostatic  and  bistatic  radar  amplitude 
patterns  are  not  the  same.  This  difference  can  be  modeled  as 

ab{x,y -u,oj)- am{x,y -u,uj)  [l  A  e{x,y  -  u,lo)\  ,  (10) 

where  e{x,y  —  u,uj)  is  an  unknown  pattern  which  contains  a  phase  function  as  well  as  a 
gain  (magnitude)  function.  Using  the  expression  for  the  bistatic  amplitude  pattern  in  (10) 
and  the  Fourier  properties  of  the  SAR’s  AM-PM  signal  in  [2,  Appendix],  the  bistatic  SAR 
reconstruction  in  the  spatial  frequency  domain  (see  (8))  becomes 

Fb{kx,ky)  =  Amiku,oj)  [l  +  E{k^„u)\  exp(jkxXn  A  jkyyn),  (11) 

^  This  phase  compensation  approximately  removes  a  shift  in  the  bistatic  data  with 
respect  to  the  monostatic  data;  the  shift  is  related  to  the  distance  of  the  two  radars  in  the 
along  track  domain  [1].  Any  residual  error  in  this  processing  also  corresponds  to  a  shift 
operation  which  can  be  viewed  as  a  component  of  the  unknown  impulse  response  h{x,y). 
Thus,  this  additional  shift  error  is  already  incorporated  in  the  general  model  of  (13);  i.e., 
h{x,y)  contains  a  smearing  component  due  to  the  calibration  gain  and  phase  errors,  and 
a  shift  component  due  to  the  residual  shift  error. 
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where 

E[2ksm6n{u),co]  =  e(xn,yn  -  u,u}),  (12) 

and  (kxiky)  are  related  to  (ku,u>)  via  (6). 

Comparing  (5)  and  (11),  we  can  write  the  following: 

Fb(^kxjky^  —  Hijtxyky^  Fmikxiky^ 
fbix,y)  ^  fm(x,y)  **h(x,y) 

where 

H{kx,ky)  =  1  +  E{ku-,i^),  (14) 

represents  the  transfer  function  of  a  linear  shift-invariant  system;  h(x,y),  the  impulse 
response,  is  the  two-dimensional  inverse  Fourier  transform  of  H{kx,ky),  and  **  denotes 
two-dimensional  convolution  in  the  spatial  domain. 

Equations  (13)  and  (14)  indicate  that  the  bistatic  SAR  image  is  a  filtered  (smeared 
and/or  shifted)  version  of  the  monostatic  SAR  image;  the  transfer  function  of  the  filter 
H(kx,ky)  is  related  to  the  phase  and  gain  differences  of  the  two  uncalibrated  radars.  This 
can  result  in  gain  and  phase  differences  between  the  monostatic  and  bistatic  SAR  images 
of  stationary  targets.  In  this  case,  the  difference  image  is  not  a  reliable  statistic  for  moving 
target  detection. 

Figure  Id  shows  the  bistatic  SAR  image  for  a  case  involving  uncalibrated  radars.  For 
this  example,  we  use  35  nonzero  coefficients  (Nx  =  5  pixels  by  Ny  =  7  pixels)  for  the 
filter  h(x,y).  The  complex  coefiicients  of  1  —  h{x,  y)  are  generated  with  a  random  number 
generator  with  a  gain  of  approximately  .2  and  a  uniform  random  phase.  Figure  le  is 
the  difference  image.  Note  that  the  difference  image  for  this  case  contciins  a  significant 
contribution  from  the  stationary  targets  (foliage)  as  well  as  the  moving  targets. 

4.  Signal  Subspace  Processing  for  Uncalibrated  Radars 
4a.  Algorithm 

The  two-dimensional  system  model  and  the  resultant  effect  which  are  described  for 
a  monopulse  SAR  system  in  the  previous  section,  are  similar  to  the  clutter  cancellation 
problem  in  one-dimensional  active/passive  imaging  problems  with  monopulse  radar  and 
shipboard  infrared  search  [3]-[6].  For  these  applications,  the  use  of  adaptive  filtering  meth¬ 
ods  has  been  suggested  [7].  One  may  consider  a  two-dimensional  version  of  such  algorithms 
for  the  problem  of  uncalibrated  radars  in  the  two-dimensional  monopulse  SAR  systems. 
However,  the  implementation  of  the  adaptive  filtering  approach  requires  inversion  of  a 


(13) 
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very  large  covariance  matrix  for  each  subpatch  of  the  target  scene,  or  a  recursive  (itera¬ 
tive)  gradient  descent  solution  [6]- [7].  Another  version  of  the  algorithm  requires  a  priori 
knowledge.of  certain  clutter  parameters  which  should  remain  unchanged  in  the  SAR  scene 
[6].  Unfortunately,  this  is  not  an  appropriate  assumption  in  a  high-resolution  SAR  image. 

Equations  (13)-(14)  indicate  that  while  fb{x,  y)  is  a  smeared  version  of  fm(x,  y),  how¬ 
ever,  the  bistatic  SAR  image  belongs  to  the  linear  signal  subspace  of  the  monostatic  SAR 
image  and  its  shifted  versions  in  the  case  of  stationary  targets.  Thus,  one  can  first  obtain 
the  projection  of  fb{x,y)  into  this  signal  subspace,  call  it  fb(x,y).  Then,  the  difference  of 
this  signal  subspace  projection  fbix,y)  and  the  bistatic  SAR  image  fb{x,y)  should  provide 
a  reliable  statistic  for  moving  target  detection. 

Let  ipeix,  y),  £  =  0, 1, 2, . . . ,  A'  -  1,  be  the  orthogonal  basis  functions  which  represent 
the  above-mentioned  signal  subspace 

^  =  0,1,2,. ..,A-1]. 

To  generate  this  signal  subspace,  one  can  use  Gram-Schmidt,  modified  Gram-Schmidt, 
Householder  or  Givens  orthogonalization  procedure  [8]-[10].  The  choice  depends  on  the 
type  of  digital  signal  processing  software  package  which  the  user  utilizes,  and  the  accuracy 
of  a  given  orthogonalization  method  in  that  software  package. 

The  size  of  the  subspace,  i.e.,  N,  depends  on  the  user’s  a  priori  knowledge  of  the 
number  of  the  nonzero  coefficients  in  the  discrete  model  of  the  impulse  response  h{x,y). 
For  instance,  if  the  discrete  h{x,y)  contains  {Nx,Ny)  non-zero  pixels,  then  we  should 
select  N  =  NxNy.  (A  similar  assignment/model  for  h{x.,  y)  is  used  in  the  adaptive  filtering 
methods  [6]- [7].)  In  practice,  the  exact  value  of  NxNy  is  not  known.  In  this  case,  an 
estimate  should  be  used  based  on  the  degree  of  gain  and  phase  differences  (calibration 
erros)  between  the  two  radars.  We  should  emphasize  that  the  algorithm  is  not  sensitive 
to  the  fact  that  N  might  be  chosen  to  be  different  from  NxNy  (as  long  as  the  difference  is 
not  very  large).  For  instance,  in  the  example  which  we  will  examine,  NxNy  =  5  x  7  =  35, 
while  we  use  A  =  9  for  the  number  of  the  basis  functions. 

The  projection  of  the  bistatic  image  fb{x,  y)  into  the  basis  function  ?/>^(a:,  y),  which  is 
identified  by  the  series  coefficient  b£  (€  =  0, 1, 2, . . . ,  A  —  1),  is  found  via  the  following: 


bt=  <  fb  ,  > 

=  h{.x,y)  ^£{x,y)  dxdy 

J  X  J  y 


(15) 


where  i/’*  is  the  complex  conjugate  of  V’-  (For  a  discrete  reconstructed  image,  the  above 
double  integral  in  the  (x,  y)  domain  is  replaced  with  a  double  discrete  sum  over  the  available 
(x,  y)  pixels.) 
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The  projection  of  the  bistatic  image  into  the  signal  subspace  is  defined  via 

N-l 

fb{^iy)=Yl  (1^) 

^=0 

The  signal  subspace  difference  image,  i.e.,  the  statistic  for  moving  target  detection  is 
constructed  via 

fdix,y)  =  fbix,y)  -  fi,{x,y).  (17) 

Consider  the  monopulse  SAR  problem  which  was  cited  earlier.  We  use  a  signal  sub¬ 
space  with  size  N  =  9  which  is  constructed  from  fm{x,y)  and  its  eight  adjacent  shifted 
versions.  Figure  If  is  the  subspace  difference  image  fd{x,y)  when  the  radars  are  fully 
calibrated.  Note  that  this  image  resembles  the  difference  image  of  the  calibrated  radars  in 
Figure  Ic.  Now,  consider  the  uncalibrated  bistatic  radar  data  of  Figure  Id  which  corre¬ 
sponds  to  NxNy  =  35  (Nx  =  5  and  Ny  =  7)  nonzero  coefficients  for  the  impulse  response 
h{x,  y).  (Note  that  the  size  of  the  subspace,  i.e.,  N’  =  9,  is  less  than  the  number  of  nonzero 
coefficients  of  h(x,y).)  Figure  Ig  shows  the  subspace  difference  image  /d(r,y)  for  the 
uncalibrated  radars.  Note  that  this  shows  a  significant  improvement  over  the  difference 
image  of  uncalibrated  radar  data  in  Figure  le.  When  we  increase  the  signal  subspace  to 
35  (i.e.,  fm.{x,y)  and  its  34  adjacent  shift  versions),  we  obtain  a  subspace  difference  im¬ 
age  for  uncalibrated  radars  which  is  almost  identical  to  the  subspace  difference  image  for 
calibrated  radars  in  Figure  If. 

4b.  Implementation  and  Computational  Efficiency 

Provided  that  the  calibration  error  function  e{x,y  —  u,uj)  is  inavriant  of  the  target’s 
coordinates,  then  the  subspace  processing  can  be  applied  in  one  step  to  the  entire  SAR 
scene.  The  same  is  impractical  for  the  adaptive  filtering  approach  [6]- [7]  due  to  the  enor¬ 
mously  large  size  of  the  covariance  matrix  which  has  to  be  inverted.  (The  size  of  the  SAR 
scene  which  is  considered  in  Figures  la-g,  is  30  pixels  by  386  pixels  which  is  a  subpatch 
of  a  larger  256  pixels  by  1024  pixels  SAR  image.  In  this  case,  it  is  even  impractical  to 
implement  the  adaptive  filtering  approach  for  the  30  pixels  by  386  pixels  SAR  subpatch 
on  a  regular  computer.)  When  the  error  function  cannot  be  modeled  to  be  invariant  of  the 
target’s  coordinates,  then  the  SAR  image  has  to  be  divided  into  subpatches  over  which 
the  error  function  does  not  vary  significantly  (which  implies  that  h{x,y)  approximately 
remains  the  same  in  that  subpatch.)  The  subspace  algorithm  can  then  be  applied  to  each 
subpatch. 

We  have  implemented  the  subspace  algorithm  using  Matlab  on  a  time-share  computer 
at  our  school.  For  the  example  of  Figures  la-g,  the  CPU  time  for  creating  the  signal 
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subspace  ^  is  5.1  sec;  the  CPU  time  to  construct  6^,  £  =  0, 1, 2, . . . ,  iV  —  1,  and  fb(x,y) 
from  (15)-(16)  is  1  sec. 

4c.  Stability  and  Noise  Effect 

Both  the  adaptive  filtering  method  in  [6]  and  the  signal  subspace  projection  of  (15)- 
(16)  seek  the  same  minimum  error  energy  solution  for  the  bistatic  SAR  reconstruction 
fb{x,  y)  in  the  linear  subspace  of  the  monostatic  SAR  reconstruction  /m(a:,  y)  and  its  shifted 
versions.  In  the  adaptive  filtering  approach,  a  solution  for  the  impulse  response  h(x,  y) 
from  the  knowledge  of  fm{x,y)  and  fb{x,y),  ^  call  it  h(x,y),  is  obtained  via  minimizing 
the  following  error  function: 

/  /  [fbix,y)-fm(x,y)  **h{x,y)Y  dxdy.  (18) 

J X  Jy 

For  the  available  discrete  data,  h(x,y)  is  modeled  by,  e.g.,  N  nonzero  coefficients  {N^ 
pixels  by  Ny  pixels,  where  N  =  NxNy).  The  resultant  solution  is  then  used  to  estimate 
fb{x,y)  via  the  following: 


fbi^^y)  =  fmix,y)  **h{x,y).  (19) 

In  the  proposed  signal  subspace  method,  the  object  is  not  to  solve  the  deconvolution 
problem  (i.e.,  to  estimate  the  impulse  response  h{x,  y))  which  is  computationally  intensive, 
and  may  yield  unstable  results  (depending  on  the  condition  of  the  covariance  matrix).  The 
signal  subspace  approach  provides  a  direct  solution  for  fb(x,y)  without  the  need  to  solve 
the  deconvolution  problem.  Furthermore,  to  reduce  additive  noise  effect,  the  user  can 
discard  the  basis  functions  which  possess  relatively  weak  eigenvalues.  (This  step  can  be 
easily  incorporated  during  the  generation  of  the  basis  functions.) 

Figure  2a-g  show  the  same  signals  as  in  Figures  la-g  when  the  data  are  corrupted 
with  a  complex  normal  additive  noise;  the  variance  of  the  noise  is  chosen  such  that  the 
signal-to-noise  ratio  (SNR)  in  Figures  2a,  2b  and  2d  is  4.5  dB;  the  additive  white  noise  in 
these  three  images  are  uncorrelated  (generated  separately).  (None  of  the  basis  functions 
are  discarded  for  this  example.)  The  moving  targets’  SAR  signatures  are  prominent  in  the 
subspace  difference  images  for  calibrated  (Figure  2f)  and  uncalibrated  radars  (Figure  2g) 
as  well  as  the  difference  image  for  calibrated  radars  (Figure  2c).  The  difference  image  of 
uncalibrated  radars  (Figure  2e)  shows  the  worst  performance  for  moving  target  detection. 

^  This  operation  is  called  deconvolution. 
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Note  that  the  K-L  transform  of  the  white  noise’s  autocorrelation,  i.e.,  the  delta  func¬ 
tion,  has  a  uniform  eigenvalue  distribution  in  and  its  complement  [11].  Thus,  the 
portion  of  the  white  noise  in  the  measurement  of  fb{x,y)  which  falls  in  ^  is  filtered  in  the 
signal  subspace  difference  processing.  However,  if  a  moving  target’s  SAR  signature  does 
not  belong  to  the  subspace  of  the  clutter  (i.e.,  they  are  orthogonal),  then  the  subspace 
difference  processing  will  not  attenuate  the  moving  target’s  SAR  signature. 


5.  Future  Plans 

We  plan  to  continue  our  theoretical  investigation  on  developing  more  compact  signal 
subspaces  of  the  monostatic  SAR  image  fm{x,y)  and  its  shifted  versions  for  representing 
the  bistatic  image  fb{x,y).  This  depends  on  the  nature  of  the  error  function  E{ku-,uj) 
which  are  encountered  in  practice.  For  instance,  in  some  radar  systems,  the  error  function 
may  take  on  the  following  form: 

E(ku,u})  =10  Eo(ku)  exp\juj  ^o(ku)]] 

the  above  two-dimensional  signal  represents  aspect  angle-dependent  phase  and  gain  errors 
that  linearly  vary  with  the  radar  frequency.  We  also  plan  to  investigate  the  implications  of 
these  error  functions  in  estimating  the  motion  parameters  of  a  moving  target  in  the  SAR 
scene. 

Moreover,  our  simulated  results  on  signal  subspace  processing  show  sufficient  promise 
that  we  believe  the  algorithm  can  now  be  tested  on  realistic  monopulse  SAR  data.  This 
study  would  reveal  the  nature  of  realistic  calibration  errors  which  could  guide  us  in  devel¬ 
oping  more  compact  and  efficient  signal  subpace  processing  algorithms. 
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Figure  lb.  Bistatic  SAR  reconstruction  (calibrated  radars) 
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Difference  Image  (calibrated  radars) 

T  1  I  I  I  I  r 


522 


524 


526 


_ I _ i_ - 1 — 

528  530  532 

Range,  meters 


534 


536 


Figure  Ic.  Difference  image  (calibrated  radars). 


Cross-range,  meters 


Bistatic  SAR  Reconstruction  (uncalibrated  radars) 
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Figure  Id.  Bistatic  SAR  reconstruction  (uncalibrated  radars). 


Difference  Image  (uncalibrated  radars) 
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Figure  le.  Difference  image  (uncal ibrated  radars) 


Subspace  Difference  Image  (calibrated  radars) 


Figure  If.  Subspace  difference  image  (calibrated  radars) 


